#!/usr/bin/env python

import numpy as np, matplotlib.pyplot as plt

rad = 4.0 # radius of the first hex ring in meters

nr = 19

th = 2*np.pi*np.arange(nr)/nr

# inner ring
x = np.cos(th)# + np.pi/6.)
y = np.sin(th)# + np.pi/6.)

plt.clf()
plt.plot(x, y, 'ro')
plt.axis([-3,3, -3,3])

np.savetxt("geometry/annulus_"+str(nr)+".txt", np.transpose((x,y)), 
           fmt='%6.3f')
